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We study stabilization of an unstable cosmic string associated with spontaneously broken 
U(1)r symmetry, which otherwise causes a dangerous roll-over process. We demonstrate that 
^ ! in a gauge mediation model, messengers can receive enough corrections from the thermal 
d • plasma of the supersymmetric standard model particles to stabilize the unstable modes of the 
string. 



1 Introduction 



The string landscape is a fascinating idea to reveal the nature of the Universe and suggests us 
that the vacuum structure of a field theory itself may be complicated. In phenomenological 
model building, this idea gives us a good opportunity to revisit supersymmetry (SUSY) 
breaking and its phenomenological applications [IrS] (See [2HZ] for reviews and references 
therein). In particular, the idea that we live in a metastable SUSY-breaking vacuum is now 
one of the most promising scenarios from both phenomenological and theoretical viewpoints. 
As is emphasized recently in [S] , if such a landscape of vacua is realized in nature, the existence 
of a certain type of solitonic objects can have an important meaning. Such a soliton can be 
viewed as an energetic impurity which causes semiclassical vacuum decay via rolling-over the 
potential hill toward a lower energy vacuum. This idea was proposed almost thirty years ago 
by several authors [SHTTj. Here we revisit it in light of the landscape of vacua and study how 
to avoid the roll-over problem in realistic SUSY-breaking models. 

In this paper, as an illustration, we consider a gauge mediation model with spontaneously 
broken U(1)r symmetry and supersymmetry. So far, various attempts have been done to build 
several types of models of spontaneously broken R-symmetry [T2H25] . One of the lessons worth 
emphasizing is a connection between metastability and large gaugino masses. To generate 
large gaugino masses in gauge mediation models, as is nicely formulated in [26] (see [27] 
for a review), one needs a tachyonic direction that leads to the SUSY vacuum in the pseudo- 
moduli space if the low energy effective theory is approximated by a renormalizable generalized 
O'Raifeartiagh model. In the models with spontaneously broken R-symmetry, the tachyonic 
direction exists at the origin [T3]. Since the U(1)r symmetry is spontaneously broken at the 
SUSY-breaking vacuum, a global cosmic string can be formed by the Kibble-Zurek mechanism 



[28l[29] at some time in the cosmic history 3 In the core of usual strings, the energy is large 
and the symmetry is restored. However, as mentioned above, at a symmetry restoring point 
in this model, there exists a tachyonic direction along the messenger direction. Therefore, the 
core of the string in this case slides down to a lower vacuum and the string transforms into the 
metastable tube-like soliton which we call R-tube. In [8], we investigated the various aspects 
of the R-tube and found that there can be light unstable modes. In these reasons, a gauge 
mediation model with spontaneously broken R-symmetry is an ideal example to demonstrate 
implications of such complicated vacua in a realistic situation (see [33] for relevant earlier 
works) . 



1 In [3D] solitons in metastable SUSY breaking vacuum were investigated and later, such solitons were used 
for cosmological constraints on gauge mediation models [31]. In addition, a metastable system with U(l) 
symmetry would decay through Q-balls. That leads to another cosmological constraint [3T| . 
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In this paper, we claim that such unstable modes of cosmic R-string/R-tube can be sta- 
bilized by the thermal potential generated by the thermal plasma of the supersymmetric 
standard model particles even with sufficiently small reheating temperature. One may think 
that as the universe expands, the temperature decreases and eventually the thermal potential 
becomes too small to stabilize the unstable modes. However, since R-symmetry is explicitly 
broken by the gravity effect, axionic domain walls connecting the R-strings are formed when 
the Hubble parameter becomes comparable to the R-axion mass and finally all string network 
disappears. Therefore, if the thermal protection of the unstable modes is valid until the time 
of their decay, it is plausible to conclude that such models are free from the disastrous roll-over 
problem. 

The paper is organized as follows. In section 2, we set up a gauge mediation model with 
spontaneously broken U(1)r symmetry. We show similarity between global string solution 
in [8] and that of the present model. We review some aspects of the string solution and 
add some new comments on the solution. In section 3, we demonstrate that the thermal 
potential generated by the thermal plasma of standard model particles lifts the potential of 
messengers and stabilized the unstable modes of R-string/R-tube. In section 4, we consider 
the model in the inflaton oscillation dominated era of the expanding universe and study the 
vacuum selection^! and stabilization of strings. In section 5, we comment on relevant issue 
and discussions. In appendix A and B, we show some technical details on R-string solution 
and tachyonic mass around it. 

2 Model with spontaneously broken R-symmetry 
2.1 Setup of Model 

Let us promote the original model studied in [8] to more realistic gauge mediation model. 
We introduce the superfields, X, <p and <fi. A pair of <fi and <fi correspond to the messenger 
fields, which have the vector-like representation, R © R, under the standard gauge group, 
SU(3) x SU(2) x i7(l). For simple illustration, here we restrict ourselves to the {7(1) gauge 
part of SU(3) x £77(2) x {7(1) as the visible sector. Then, we consider the messenger fields, 
<p and <p, which are singlets and have the {7(1) charges 1 and —1. Its extension to the full 
£{7(3) x £{7(2) x {7(1) gauge theory is straightforward as we will give a comment later. 

We assign R-charges to these superfields as, R[X] = 2 and R[<f>] + R[4>] = 0, and the 
superpotential is obtained by 

W = KXcpfi- pt 2 X. (2.1) 
2 See }34f(58] for early attempts for vacuum selection by exploiting the thermal potential. 
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Note that we introduced k to control the vacuum selection. As is used in J39JII0] by taking 
k to be small, the tachyonic mass \m^\ along 0, can be small, which is favorable for the 
realistic vacuum selection. We consider the effective Kahler metric, 
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which are similar to those in [8]. We take all of parameters, k, /i 2 , m 2 and A, to be real and 
positive. This model possesses the SUSY vacuum with the moduli space, 
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and the SUSY breaking vacuum, 
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where the U(1)r symmetry is also broken. The former is the true vacuum, while the latter is 
the metastable vacuum whose vacuum energy is V = \ft (l — 

For later convenience, let us introduce dimensionless variables by 
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then the Lagrangian is of the form 
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where the ellipsis denotes the fermionic partners and the supersymmetric standard model 
particles including gauge fields and matter fields, and V(A) is given by 



\A 
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(2.7) 



Basic ideas we will show below can be simply demonstrated by this simplified model. Since 
throughout this paper, we do not consider a soliton with the winding number of U(l), we can 
take the vacuum expectation value (VEV) of gauge fields vanishing in constructing solitons. 

Various quantities appeared in the first line of (12. 6p are characterized by two dimensionless 
parameters A and e. For instance, the existence of the SUSY breaking vacuum requires 



1> Tx >0 - 
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e 2 4A 



(2.8) 
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Later, we consider the vacuum selection in the early stage of the Universe. Since in the 
early Universe field values are assumed to be around the origin, if the tachyonic mass of X, 
\mx\, is larger than one of the messenger field, \m^\, one may expect that the supersymmetry 
breaking vacuum is preferable. In the dimensionful representation, these two tachyonic masses 
are given by 

,2 h 1 2 ,.,,2 



X = ~2m^' < = - K ^- ( 2 - 9 ) 



Thus, the following inequality 



m\ 2 1 
— > I ±± _ < _ 

ml e 2 2A 



> 1 — < —r, (2.10) 



is required for selecting the SUSY breaking vacuum. Note that, according to the reparametriza- 
tion (12. 5p . a small k with fixing m, p corresponds to a large e. As discussed later in section HJ 
the thermal effect turns out to loosen this condition. 

For later convenience, we show the gravitino and axion masses, 



mv^VrW^-, m a = ^m 3/2 (g-\X ))-^^, (2.11) 

where g -1 (Xo) = 1 — 4- and M p \ denotes the reduced Planck mass. Using the gravitino mass 
m3/2 with dimensionless parameters A, e, k, many scales can be rewritten as, for instance, 
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48A 3 /V12A - 3 m 3/2 18 2 (4A) 3 / 2 V / 4A^T m 3/2 

_ x — — ~ 2.1 x 10 x — - — x — — . (2.12) 

m\) e 2 n M p i e 2 k GeV V ; 

In terms of (12. 8p and the scale of the gravitino mass in gauge mediation, we observe \mx\ 3> 
m a in wide range of parameter space (A, k). Therefore, we assume the condition in this paper. 

2.2 R-string review 

Here we quickly review the R-string studied in [8]. Since the R-string solution satisfies the 
D-flatness condition, the equation of motion for R-string becomes exactly the same as the one 
studied in [S]. Hence, we simply see some of aspects. 

The R-string solution without a hole inside, s = s = 0, is a solution of our model. The 
R-string solution A(x^) = A^ 1 (p)e ind , (pe md — X\ + ix 2 ,A s ° l E K>o) with a given winding 
number n is defined by imposing the equation of motion derived from the following 'reduced' 
action, 

S = ~2nl dpp{^((^-) 2 + ^)+V(A)l (2.13) 
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with boundary conditions A(0) = and A(oo) = 1. The R-string solution corresponds to the 
solution of the system described by the above Lagrangian and we can find it numerically (see 
appendix |X] for details of the solution.). However, as shown in [5] such an R-string is unstable 
and transforms into an R-tube with non-zero s = s inside. Since our main motivation is to 
stabilize such unstable modes, let us quickly review the analysis to get "masses" of the modes. 
To see that, let us consider a linearized equation for s around the R-string solution, 

" ~ (p^) + ^pot(^ ol (p)) a = q 2 n s, V pot (A) = -|V(A) + ^\A\ 2 . (2.14) 

Here, the eigenvalue q\ for the R-string solution with the winding number n depends on A 
and e, i.e. q\ = q 2 (X,e), and, an observation < A h r f < 1 tells us the lower bounds of V pot 
and q\ as 

2 2 2 

Taking the limits of A — > oo or e — >■ 0, we can show that 



V pot (A)>V pot (0) = --, ql>— 2 . (2.15) 



lime 2 ^ = -2, lim^ = -2, lime^ = -2. (2.16) 

n— >-oo A— s-oo e— >u 

For the proof of these facts, see AppendixjBl Furthermore, we experientially find that, with 
several typical sets of parameters A and e, 

< -e 2 q\ < -e 2 q 2 2 < ■■ ■ < 2, (2.17) 

and interestingly 

- e 2 gi 2 < 2 (l - <-e 2 ql (2.18) 

although we have no proof. We numerically solved this equation for several parameters as 
shown in Figure [1] and [2j The left and right panels in both figures correspond to n — 1 and 
2, respectively. 

For later convenience, we introduce the "physical" mass for the unstable mode. The 
eigenvalue q\ for the R-string with winding number n is translated in terms of the mass for 
the canonically normalized messenger field as 

< = q l(^f) =^xV- (2.19) 

We will use it later for the physical considerations. Here an non-trivial fact m 2 ^ = = — nfi 2 
is observed. 
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Figure 1: Tachyonic masses for unstable modes of n = 1 and n = 2 R-strings with X = 0.27. 
Solid lines represent e 2 q 2 = —2(1 — l/v4A) and a dashed line in the left panel is e 2 q 2 = 
—2 x 6 x 1CT 7 and dashed lines in the right panel are e 2 q 2 = —2, —2 (l — 0.47e 2 ^ 3 ) . 




Figure 2: Tachyonic masses for unstable modes of n = 1 and n = 2 R-strings with A = 1/2. 
Dashed lines in the left panel are e 2 q 2 = —2, —2 x 0.065 and dashed lines in the right panel 
are e 2 q 2 = -2, -2 (l - 0.39e 2/3 ) . 
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3 R-string in thermal plasma 



Now we turn to the stabilization of the R-strings. As shown in the previous section, R-strings 
in the (Minkowski background) vacuum are unstable and unavoidably deform to R-tubes. As 
discussed in [8], R-tubes are basically unstable objects and cause, in turn, the roll-over process 
to push the whole Universe into the unwanted SUSY vacuum with a finite time. Therefore, 
one may consider that the SUSY-breaking models with the spontaneous R-symmetry breaking 
are disfavored if the R-symmetry has once been restored in the cosmic history and R-strings 
have been formed. 

However, the realistic Universe is neither described by the Minkowski space time nor in 
the vacuum. It contains several ingredients that can deform the potential for the R-string 
sector. It may be possible to stabilize the R-strings. The most effective contribution would 
be the thermal potential in the existence of the thermal plasma and hence we focus on the 
thermal effect in the following §. 

Here we should note that we do not have to stabilize R-strings up to the present time. As 
explained in [H] , R-string networks, if they are stabilized, vanish when the Hubble parameter 
of the Universe becomes sufficiently small, H ~ H a = m a [32] • This is because the constant 
term is added in the supergravity superpotential to compensate the vanishing cosmological 
constant, and such a constant term explicitly breaks the R-symmetry and makes the R-string 
healthily unstable. From the explicit R-symmetry breaking effect, the R-axion, which is the 
phase direction of the SUSY-breaking field acquires a mass. Thus, when the Hubble parameter 
decreases sufficiently, the R-string network turns to the R-string-domain wall system, which 
collapses to R-axion particles immediately [32]- Therefore, we only have to stabilize R-strings 
up to the Hubble time corresponding to H = H a . 

We also note that high reheating temperature is not necessarily required in our discussion. 
In the gauge mediated SUSY-breaking, there is a severe constraint on the reheating tempera- 
ture from the gravitino proble nfl Therefore, thermal plasma seems to be difficult to modify 
the potential for the R-string sector. However, even before reheating during the inflaton oscil- 
lation dominated era, thermal plasma of supersymmetric standard model particles does exist 
from the partial decay of inflaton quanta, of which temperature is larger than the reheating 
temperature. Thus, it is plausible to consider the thermal potential for the R-string sector 
that is generated by the thermal plasma. In this section, we demonstrate that the unstable 



3 Thc Hubble induced masses generated from the Planck suppressed interaction between the R-string sector 
and inflaton also exist during the inflaton oscillation dominated era before reheating. However, we can show 
easily that it cannot stabilize R-strings sufficiently in the wide range of parameter space. 

4 For the cosmological constraint from R-axions, see [41] . 
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modes shown in the previous section can be stabilized by the thermal effect assuming the 
existence of thermal plasma with sufficiently high temperature. We will revisit the cosmic 
history with R-strings in the next section. 



3.1 Thermal potential for messenger 

In this subsection, we study the thermal potential in the hidden sector generated by the 
supersymmetric standard model thermal plasma through the gauge interaction. Here we 
assume that the temperature of thermal plasma is relatively low, T < kXq and the moduli 
field and messenger field are not thermalizec0 since we are interested in the last stage when 
the stabilization mechanism is effective. 

Since nonvanishing messenger field values generate the Standard model gauge boson (and 
gaugino) mass, their one-loop effective thermal potential generates thermal correction to the 
scalar potential for the messenger fields (see (13] for a review and references therein), 

Vfcennaito 0) - 3^ J s (m|/T 2 ), m\ = ^ (|0| 2 + |0| 2 ) , (3.1) 
where factor 3 comes from the degree of freedom of the massive U(l) vector bosons^ and 



POO 

J B ( y ) = dxx 2 log 1 - e"V^ . (3.2) 







For smaller messenger field values, we can use the high temperature expansion, 

4 2 2 / 2 \ 3/2 -i 4 2 

2 /rp2\ 71 71 m B 71 ( m B \ 1 m B m B 



Mm%/T°) = - log 



45 12 T 2 6 \ T 2 J 32 T 4 ° a B T 2 

— g(-^(- + 8(^r. <-) 

with loga^e = 5.4076. Then, the messenger fields acquire the effective mass, mT, as 

Khe rm ai(0,0)^f^T 2 (|0| 2 + |0| 2 ), =► mJ = |T. (3.4) 

We can expect that the R-string is stabilized if the thermal mass is larger than the physical 
mass of the tachyonic mode discussed in the previous section. The ratio of the messenger 



5 For smaller temperature, the number density of heavy fields is Boltzmann suppressed, if ever, and hence 
they cannot be in thermal equilibrium in the expanding Universe. The mass of messenger fields around the 
R-string core is typically small enough to be thermalized. However, the width of R-string is too small for the 
messenger fields to be thermalized and we do not consider the effect of messenger thermalization. 

6 Rcmcmbcr that for simplicity we consider only the singlet messenger fields 4> and 4> with the U(l) charges 
±1. It is straightforward to proceed the argument to the full SU(3) x ST/ (2) x U(l) theory. 
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mass to the tachyonic mode mass is given by 



\m n \ e 2 \q n \ 



(3.5) 



which is larger than one in a wide range of parameter space. Thus, even in a low temperature, 
T < kX , the thermal mass can be comparable to the physical tachyonic mode mass for the 
R-string, and hence it is enough to use the above thermal potential to study stability of 
tachyonic modes. 

To complete the argument on thermal effects, we comment on the higher temperature situ- 
ation, T > kXq. In this case, the messenger fields are also thermalized and extra contribution 
to thermal effective potential to messenger direction is generated, 

V^4^(J B (m 2 B /T 2 )-J F (m 2 B /T 2 )), m| = £ (|0| 2 + |0| 2 ) (3.6) 



poo 

\ dx x 2 log 


1 + e ~^ x2+v 


/o 





J F ( y ) = dxx 2 log 1 + e-V^+w . (3.7) 



where 



Here, the fermion contribution comes from the diagram with gaugino and messengino, whereas 
the additional scalar contribution comes from the messenger scalar loop, which acquires mass 
from the D-term. We should note that the thermal potential for X is also generated by 
messenger fields, V ~ k 2 T 2 \X\ 2 + However, if k is small, the thermal potential can 
be much smaller than the zero-temperature mass of X at T > kXq. Thus, we can neglect 
thermal corrections to the potential for the X field here. 

Finally, we comment on the extension of the gauge theory. So far, we have considered 
the singlet messenger fields, <fi and 0, with the U{1) charges ±1 for simplicity. It is straight- 
forward to extend our discussions for generic messenger fields, <fi and <fi with the vector-like 
representations, R© R, under the SU(3) x SU{2) x U{\) gauge symmetry. In such a generic 
case, all of the gauge bosons contribute to induce the thermal potential of the messenger 
fields. Such a thermal potential is obtained by replacing g 2 in (13. ip and ( 13 .4p by 

glff = s 3 2 cf + glc?\K) + g 2 q % (3.8) 

(3) (2) 

where C 2 (R) and C 2 (R) denote the quadratic Casimir indices of the R representations of 
SU(3) and SU(2), respectively, and denotes the hypercharge of 0. Thus, our discussions 
in the following sections can be extended into the full visible SU(3) x SU(2) x U(l) gauge 
theory by replacing g 2 — > g 2 ^. 
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3.2 Tachyonic mode around R-string 

Now we are ready to study thermal effects on R-string. Firstly, let us see how the tachyonic 
mass of the R-string can be modified by the thermal effect in (13.11) . To study this, we have 
only to pay attention to an infinitesimal fluctuation |0| 2 = \<p\ 2 = p 2 \s\ 2 / k <C p 2 /n, which 
indicates that we can use the high temperature approximation, 



30 p 1 + 8^p 2 



thermal = p" + f o |*f + " " " • (3-9) 



The second term in the above gives a constant shift of the potential in (I2.14p as 

g>_T 2 
8ne 2 /i 

and resultantly a mass eigenvalue q 2 (T) also accepts the constant shift 



V pot (A) -+ y P ot(A,T) = V pot (A) + ±^— (3-10) 



q 2 n (T) = q 2 n + -f^ 1 ^. (3.11; 

It is straightforward to define the n-th critical temperature, T n , where q 2 {T n ) = 0, that is, 



4 / e 2 q 2 4 

T n = -xd ^ x up 2 = ~m n , (3.12) 

9 V 2 g 

and then we find the following properties 

< Ti < T 2 < ■ ■ ■ < -^^p = T^, lim T n = -^^p, limT n> i = -^^p. (3.13) 

g A^oo g e^O g 

At a temperature T larger than T n , the tachyonic mass changes to massive one q 2 {T) > and 
therefore the R-string with a winding number n is stable. 

As an illustration, we show the low energy effective potential for the light mode. As in [8] , 
to uncover the existence of light unstable mode, which is almost frozen in the relaxation 
method, we vary the initial profile function. To be concrete, the following initial conditions 
with various values of po are used for A = f(p)e me , s = h(p) 

tt \ 1 +tanh(2(p-p )) 

f(P) = 2 tanh(p), (3.14) 

= 1 - tanh(2(p - po)) 1 + tanh(2(p + p )) 
[P) 2 2 

Also, to represent the size of the tube, we define the following parameter, 

j™dpp 2 s'{p) 2 

Rs = roo , ., n 2 • (3.15) 

Jo dpps'ip) 2 



10 



We also introduce the dimensionless temperature as 



where as is shown in (13. 3p . the parameter for the thermal potential is written as 



i z 4 4 

1 m% 2 _ # 4 



^ < oo., ,.u, (3-17) 



4tt 2 T 2 32vr 2 e 2 K 2 £; 2 32vr 2 e 2 fi: 2 A; 2 ' 
To see how the stabilization mechanism works, we define the dimensionless energy as 

E(r) = A. E -f^ = _ i V (l)A 2 , (3.18) 

v ; m 2 2ttAz m 2 27iAz 2 v ; ' K J 

at the relaxation time r (for quick introduction to the relaxation method, see appendix A 

in [8]). Here we removed a contribution of the vacuum energy density E vev . A denotes the 

cut-off of the energy and we set A = 50. Since we cannot follow the evolution of the system 

in the relaxation method to r — > oo, we evaluate it at r = 50. In Figure El we plot R s 

and E(t = 50) for various initial conditions for k = 0.001 and k = 0.01. Here we take the 

parameters as A = 0.27, e = 1, and n = 1. For these parameters, we can expect the R-string 

is stabilized for k > y^Ti/T^ 0.00155. As one can see from the Figure, R s = is the 

local energy minimum of the field configuration, which implies the R-string is stabilized for 

k = 0.01, whereas it does not seem to be local minimum for k = 0.001, as expected. Thus, 

we conclude that the R-strings are stabilized for T > T n . It is interesting that the R-strings 

are stabilized even if the thermal mass for the messenger field ~ gT is smaller than its zero 

temperature mass at the origin, \m^\ = y/^fi. Finally, it is worth emphasizing that in the 

Figure [31 we draw the effective potential until relatively larger value of R s just for reference 

where high temperature expansion is no longer valid. Thus, only the behavior in small R s 

region is reliable, which is enough to guarantee the stability of the R-string. 

3.3 Uplifted SUSY vacuum and stable R-tube 

We here point out that there is a phenomenologically safe situation even when the R-strings 
are unstable and they deform to R-tubes. One of the interesting features of thermal effects 
is lifting the lower energy vacuum. Since SUSY is broken by the thermal effect, the SUSY 
vacuum can be lifted and the SUSY-breaking vacuum becomes thermodynamically favored 
compared to the lower energy vacuum. To see that, let us focus on the messenger direction 
at X = and |0| = \$\. 

*W = (^ 2 - V 2 V + I hi I . (3 - 19) 
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Figure 3: The low energy effective theory for A = 0.27, e = 1 and n — 1. Left panel is for 
almost zero temperature k = 0.001. Right panel is for k = 0.01. 



where = g\<p\ = g\<j)\ and 



J(x) 



2tt' 



;(J B (x) - J B (Q)). 



(3.20) 



Here, we subtracted Jb(0) because the SUSY breaking vacuum also gets lifted by thermal 
effect. It is given by constant term at <p = <p = 0. Thus, it is useful to subtract it and define 
the potential as above. Here we can show the following inequalities 



< xj'(x) < J(x). 



The would-be SUSY vacuum is determined by 



dV 

w v mgr 

dip 2 



H 



T 2 J' = 0, 



k_ T 
9 V 



(3.21) 



< (3.22) 

K 



The critical temperature, T cr , at which the would-be SUSY vacuum and SUSY breaking 
vacuum become the same energy density, is determined by the following equation^; 



— ), => T = T cr (-,X 



H 



(3.23) 



By taking derivatives with respect to k and A of the above equation, we find that T cr (K,/g 2 , A) 
is a monotonically increasing function of k and A as 



dT cr 



4A 2 



dV r 



mgr 



dT 



If 



8T 

> u - 1 - cr 



dV r 



iu«r 



dT 



T— T C r,' 



T=T r , 



>0, 

> o, 



(3.24) 



7 Note that the field configuration in the realistic expanding Universe does not follow perfectly the R-tube 
solution. However, at least the R-tubc solution is stable, the dangerous roll-over process will not cause and 
hence this consideration gives a good condition to avoid the roll-over problem. 
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where we used the inequalities (13.21]) . 

For A > 1 or k > 1, that is 2 ac <C 47r 2 T c 2 r , the potential reduces to 



V mg ^[j^ 2 -^ + \t^\ (3.25) 



and thus, the functions vac and T cr are explicitly given as 



. -£U-€£\ = M T„ = — fix , l l —. ,3,6) 



vac 



Then, the inequality ( 12 . 1 8[) means that, for sufficiently large k, 

Ti < < T 2 . (3.27) 
In the opposite limit, A ~ 1/4 or k <§C 1, we also find the explicit formal 

K,^(^-^) 2 + ^ ^(^-i-))*. ,3,8) 

For \fRj 1 g <C 1, therefore, we find the inequality T n < T cr with an arbitrary n. 

Here, we show that the critical temperature of the R-string with the minimal winding 
number is always lower than T cr . By using (13.221) . ( 13.241) and (13.261) . we fine the following 
relation 



dT cr T cr AJk 

From the properties noted in the above we can draw a phase diagram as Figure HI Note that 
an R-tube with the minimal winding number n — 1 is always unstable. 

As an illustration, we show a stable R-tube solution with the winding number n = 2. 
We take parameters k — e — 1 and k = 0.5, which allow us to use k ^> 1 limit as ( 13.261) 
though it is not good enough. From the Figured], we see that the tachyonic mass is |m 2 | ~ 
(1.03/v^) x for A = 0.27 and e = 1. Since k = 0.5 < 1.03, the R-string with n = 2 is 
unstable and turns out to become the R-tube in that case. On the other hand, k = 0.5 is 



enough to uplift the SUSY vacuum to stabilize the R-tube as k = 0.5 > J 2 - 1/y/X = 0.275. 
In the Figure[5]we show the low energy effective potential for the R-tube solution with k = 0.5, 
A = 0.27, e — 1. Clearly there exists the minimum around R s ~ 7.2. 



8 In our model, the gauge group is completely broken in the SUSY-vacuum and hence thermal logarithmic 
potential [44) for low temperature regime would not arise. However, even if it appears, it also has an effect to 
lift up the SUSY-vacuum, and hence the basic feature docs not change if it is positive. In the case of negative 
thermal logarithmic potential, the stabilization mechanism for R-tube does not work with A ~ 1/4 or k < 1. 
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T 




Figure 4: The left panel shows phases of R-string/R-tube. With a winding number n, an 
R-string is stable even including quantum effect at high temperature T > T n ,T cr , and is 
unstable in any sense at T < T cr , T n . For T cr > T > T n , an R-string is classically stable and 
for T n > T > T cr an R-tube is stable. 



Finally, let us comment on the validity to use the thermal potential. In deriving the 
thermal potential, we assumed homogeneity at least of order 0(T~ l ). On the other hand, 
when we study the stability of solution we put the R-string/R-tube solution which generates 
inhomogenioety in the space. A condition of the validity to the thermal potential can be 
roughly estimated as 

9 4 < T ^ ^ds < T. (3.30) 

In the stabilized R-string solution, there is no space dependence of s. Thus this condition is 
satisfied. However, as for the R-tube solution shown above, there is space dependence. Using 
the dimensionless temperature in A3. 16j) . we can estimate the validity condition, 

4 

ds < -k. (3.31) 
9 

The numerical analysis on the R-tube solution for k = 0.5 leads to ds ~ 1/4. Thus, our 
numerical results are still under control for a proper value of g. 

4 Application to the expanding Universe 

Now we study how the discussion in the previous section is applied to the constraints on the 
realistic cosmology with expanding Universe. Here we focus on the vacuum selection and the 
way to avoid the roll-over problem. There are also constraints from the cosmological gravitino 
or moduli problem, and R-axion particle abundances, but these are beyond the scope of the 
present paper. For the discussion on such problems, see [ID] and [4Tj . 
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Figure 5: R-tube solution with n = 2, k = 0.5, A = 0.27, e = 1. A minimum exists around 



R s ~ 7.2. 



Let us follow cosmic histories in our scenario schematically. The thermal mass for the 
messengers is given by mT = c^gT, with being the numerical factor of the order of unity 
that counts the number of fields contributing the thermal mass. If messengers enter thermal 
plasma, it generates thermal mass to the X field, as m x = cx^T, where cx is the numerical 
factor of the order of unity. Note that during the inflaton oscillation dominated era, the 
temperature and the Hubble parameter are related by 



T 



1/4 



90 

n 2 g* 



(T 2 R M pl H) 



1/4 



(4.1; 



where g* is the effective relativistic degrees of freedom and Tr denotes the reheating tempera- 
ture. There are also the so-called Hubble induced mass for X and messengers, 



~ rrv\ ~ 



which is generated from the Planck suppressed interaction between inflaton and these fields 
in supergravity. For small enough k, 



5 

K < I - 

2 



1/4 



n 2 9* 
90 



1/8 



\m x \ 



1/4 



(4.2) 



the Hubble induced mass overwhelms the thermal mass for X fields when H > \mx\, and in 
the following, we consider such a situation. Here, we assume the inflaton oscillation dominated 
era since the gravitino problem requires relatively small reheating temperature. 

At a high temperature in the inflaton oscillation dominated era, the thermal masses or 
the Hubble induced mass for X and messenger fields restore the symmetry and these fields 
are set at the origin. As the temperature and the Hubble parameter decreases, the Hubble 
induced mass or the thermal mass can no longer fix the fields at the origin and symmetries 
are spontaneously broken. If the Hubble induced mass for the X becomes inefficient earlier 
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than the thermal mass for messenger fields, the SUSY-breaking vacuum is naturally selected 
associated with the R-string formation. This condition can be written as 

/2\ 1/4 / 90 \ 1/8 

K| = yfcn < m T JH = m x ) = ^g[-\ — {T 2 R M pl m x ) ll \ (4.3) 



which is rewritten in terms of the constraint on the reheating temperature 

T R > T A = 



5ir 2 g, x 



(4.4) 



72 J \{c^gYM pl m x 
In the other case, if \mx\ > \ m <b\, the Hubble induced mass for X always becomes inefficient 



before that for cj), become inefficient as discussed in section 12.11 Therefore, one of the 
inequality (12.10)) and the inequality f l4.4j) is needed to be satisfied for selecting the SUSY- 
breaking vacuum as showed in the Figure |6j 



r 2. 
Log,,, — 




10 4A 



Figure 6: The region satisfying one of the two inequality T| r > T4 (defined in (14.12ft ) and 
\mx\ > \ m (/>\ with (m3/ 2 /GeV)/t -1 / 3 = 10~ 6 , 10~ 3 , 10° which are colored by red, blue, and 
green respectively. The black region is excluded by the condition of the existence of the 
SUSY breaking vacuum. 



After X acquires the nonvanishing field value, the R-strings are formed and our discussion 
in the previous section can be applied. We have shown that any R-strings with arbitrary 
winding numbers are stabilized for sufficiently high temperature, 

T > T^. (4.5) 

However, practically, it may be sufficient to protect only the mode with n — 1. In this case, 
the condition for the stabilization can be written as 

T>mm(T 1 ,T cr ) = T 1 . (4.6) 
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As shown in (I3.29p . T cr is always higher than 7\. Since the R-string networks turn to the 
R-string-domain wall networks and immediately decay to R-axion particles at H = m a [42] . 
we can avoid the roll-over problem if the stabilization mechanism discussed above works at 
that time. As shown in the previous section, for y/n/g <C 1, T n < T cr is satisfied. Thus, the 
above constraints can be rewritten as follows, 



/2\ 1 ( 90 \ 1 

T(H = m a )>T 00 ^^fi<ml(H = m a )=c^gi-\ f^-J (T 2 M p im a ) 1/4 , (4.7) 
which is rewritten in terms of the constraint on the reheating temperature, 



TR>TBs[ ^y t (^J) i = MT A , (4.8, 
72 / \{c^gpM p im a J V m a 



and 



e\a I /2\ 1 ( 90 \ 1 

T{H = m a ) > T x & -j^V^lJi < m^{H = m a )) = c^g f -J f —J (T^M pl m a ) 

which, in turn, is expressed as the constraint on the reheating temperature, 



1/4 



(4.9) 



T R >T c = 



2^2 



T B . (4.10) 



1152 ; V(c^) 4 M pl m a/ 

Note that a quantity |e 2 qf/2| 4 (< 1), which is compared with a value in Eq. fl2.12p . can take a 
value 1 in the limit of A — > oo and also 10 -22 in a case of A = 0.27, e = 1, and thus there are 
both cases of Ta > T c and T4 < T c whereas inequalities T B > Ta,T c are always satisfied. 

In Figure [7J we illustrate thermal histories by showing Hubble parameter dependences of 
various mass terms. In this figure we assumed Ta > Tc, \mx\ > \ m rj>\ just to draw this figure 
concretely. The blue line is the Hubble induced mass linear in H. The red lines represent 
the thermal masses for the messengers with three different reheating temperatures. The 
horizontal green lines are the sizes of the tachyonic masses \mx\ and |m^| at the origin with 
the following relation, 

^ = ^ (4 11) 

Note that depending on the two parameters, the order of the two masses is changed. Below 
we will show that the cases 1 and 2 in Figure [7] realize the successful vacuum selection. At 
the point A, the Hubble parameter becomes comparable to the size of tachyonic mass of X 
at the origin, so the vacuum starts to slide down to the SUSY breaking vacuum. Sine the 
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condition H4.3[) is satisfied in these two cases, the SUSY breaking vacuum is successfully chosen. 
The points B and C represent the points where T(H = m a ) = and T(H = m a ) = 7\ 
are satisfied, respectively. In summary, the case 1 passes both the vacuum selection and 
stabilization of all R-strings with arbitrary winding number, whereas the case 2 passes the 
vacuum selection but stabilizes only the R-string with n = 1. The case 3 is ruled out by the 
vacuum selection if |m^,| > \mx\- With typical values of A, e, the condition f 14 . 8 [) is too strong 
to be satisfied, and the case 2 or the case 3 are preferred where a value of e 2 qf/2 plays a quite 
important role. 

log Mass 

Hinf 

\m x \ 
l™zl 

ma 

H~m x H~m a -log H 

Figure 7: In the case 1, the reheating temperature is in Tr > Tb and the SUSY breaking 
vacuum is preferable. R-strings with all winding number are stable until their decay. In 
the case 2, the reheating temperature satisfies Tb > Tr > Ta > Tc ■ The vacuum selection 
is successful but some of higher winding R-strings are unstable. In the case 3, reheating 
temperature is lower Ta > Tr > Tc ■ The vacuum selection is unsuccessful when the tachyonic 
masses at the origin are \m^\ > \mx\- 

As a reference, we compare these results to the constraint on the reheating temperature 
from the gravitino problem [35] , 

T R <T f ,6.4xlO° m3/2 y^)- 2 , (4.12) 

where m\ denotes the gaugino mass. This condition comes from the constraint for the ther- 
mally produced gravitinos at the time of reheating not to overdose the Universe. Note that 
there may be late time entropy production from the moduli decay, which would relax the 
constraint on the reheating temperature, but this upper bound on the reheating temperature 
is a good reference value. 




18 



Figure [8] shows the window satisfying T| r > T4 in the parameter space. We can see that 
there is a definite parameter space. For relatively large k, a small gravitino mass is favored 
by the vacuum selection. Let us show the existence of non-zero window for T| r > Tg. In 
Figure O we show regions satisfying T| r > T#. Here the coefficient c^g is roughly estimated 
order one, c^g = 0(1). We adopt the gaugino mass in the gauge mediation as 



m x = 0(0.1) x 



1 g- l (X )V\^ 



167T 2 



Km 



(4.13) 



where we assume a small hierarchy between gaugino mass and squark mass by introducing 
the small parameter 0(0.1) coefficient which can be interpreted as the R-breaking effect for 
Majorana gaugino mass. Finally we show examples for the case T| r > Tc in Figure [101 



e=l/5, A=|0.27 




Log,,, [m w /GeV] 



Figure 8: The regions satisfying the condition T B > T| r > T A (> T c ) that correspond to to 
the case 2 in the Figure\7\ We take g* = 220 by assuming the messengers are not thermalized. 



Here we have used the numerical results for the tachyonic masses shown in Figure [H 

In conclusion, we surely have the scenarios that are free from the roll-over problem, avoid- 
ing the thermally produced gravitino problem. For this mechanism the value of a tachyonic 
mass of a fluctuation around the R-string configuration, qf, often plays a quite important 
role. This value is very sensitive with details of the model and can take a value with a quite 
wide range as we discussed, and thus, one has to always check this when building a more 
realistic model. Of course, there are also constraints from the moduli problem, the moduli 
induced gravitino problem, and R-axion problems, which are beyond the scope of the present 
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0.0 0.2 0.4 0.6 0.8 1.0 4 A Log,„ [m ss /OoV] 



Figure 9: The regions corresponding to the case 1 in Figure\7\ T| r > Tb- We assume g* = 220. 
All winding strings are stable until the times of their decays. In the left panel the regions col- 
ored by purple, red, blue, and green correspond to (m 3 / 2 /GeV)K~ 1 ^ 5 = 10~ 9 , 10~ 6 , 10~ 3 , 10°, 
respectively. A region with \mx\ > is almost excluded by this strong condition. 




-8 -6 -4 -2 

Log ]0 [m 3;2 /GeV] 



Figure 10: Three colored regions satisfy Ta > T| r > T c corresponding to the case 3 
in Figure We took X = 0.27 and = 220. Tachyonic masses we used (qf , e) = 
(10- 5 - 4 ,1.4),(10- 4 - 8 ,2),(10- 4 - 3 ,2.4). 
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study. We emphasize that our study is helpful for the future model building and studying the 
constraint on such models. 

5 Discussion 

In this paper, we focused on gauge mediation models and studied the thermal potential to 
the messenger generated by the thermal plasma of the standard model particles. However, 
in the gravity mediations, interactions to the standard model particles are suppressed by the 
Plank mass. It is not obvious to get large enough thermal potential to stabilize the unstable 
modes. However, in this case, there is one interesting possibility to stabilize unstable modes. 
Supergravity corrections also provide a positive mass term to the fields and 0, 



where / = /i 2 is the F-term of X and c denote the constant term in the superpotential, we 
have not explicitly written similar terms for and we used the condition for the cancellation of 
the cosmological constant. From the last equality, we see that the mass is positive. In gravity 
mediation models, the SUSY breaking scale is relatively large. Therefore, if the gravitino 
mass satisfies the following inequality, 



then it may be possible to stabilize R-strings with low winding number. As long as k is very 
small, this condition is satisfied. Also, it is worth mentioning that in such a large SUSY 
breaking, the axion mass can be larger than the symmetry breaking scale, m a > \mx\- In 
this case, explicit breaking effects are significant, so the R-string is immediately broken after 
its formation. 

From our study, it may be plausible to assume that in a realistic gauge mediation model 
with spontaneous R-breaking, rolling over the potential hill by the unstable solitons does 
not occur. In this case, we have to impose further cosmological conditions, such as the 
moduli problem or the moduli induced gravitino problem studied in Ref. [30] • Moreover, as 
intensively studied in [31], when the R-string decay by axionic domain wall, it produces a 
large amount of R-axions. As for the long-lived R-axion, cosmological constraints such as Big 
Bang Nucleosynthesis and CMB observation severely constrain the parameter space [31]. 






(5.3) 



m n \ < m 3 / 2 . 



(5.4) 
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A Property of a n 

Near the core of R-string, the solution behaves as 

K° l (p) = Kp) n for a„p«l. (A.l) 

For large A a property of a coefficient a n can be discussed as follows. Let us consider the 
following limit of the action after extracting the vacuum energy 

Km § = -27T £°dpp | [{^f) 2 + ^ 2 ooj + \(l - ^o) 2 j + constant, (A.2) 

with a redefinition of the function and the coordinate 

Aoo(p) = lim A(p = V\p) . (A.3) 

A^oo \ / 

In this limit, we also have a string solution A^lp) = f n (p) with winding number n and this 
solution behaves as 

^oo(p) = fn(p) ~ {KPT for a n p < 1. (A.4) 

Note that a n is a constant and depends on only the number n. We know that the radius of 
R-string is estimated as p ~ 2n, and this fact implies that 

/ n (2n)=0(l) a n = -xO{l). (A.5) 

n 

Figure [11] shows the numerical result in the limit of A — > oo f ]A.2j) , where a n is well approxi- 
mated by a n = 1/ V4.5n + l.ln 2 . Therefore, we conclude that A dependence of a n for large A 
(and sufficiently large n) is given by 

a n = ^= = — — x C(l) for A > 1. (A.6) 
V A nv A 
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Figure 11: Dots are numerical data of a = a n and a solid line gives a = 1/y/A.bn + l.ln 2 . 



B Properties of the tachyonic mass 

The first excited mode around the R-string has the point p = p c so that 



d f ds 



dp \ dp 



P=pc 



0, & V pot (AT(Pc)) = q 2 n- 



(B.1) 



This value of p c must represent a scale of this excitation mode and thus the following quantity 
should be of order one, 



7 B =^_ Vpot (0)xp c = O(l). 



(B.2) 



The above two equations give rough estimations of g 2 , p c as follows. Let us assume that the 
following inequality 



e\l 2 1 

— ^ < 1 . 

2 e 2 2A 



Under this assumption fIB.lj) can be solved as 



(B.3) 



(B.4) 



Therefore we can use A^ ol (p) « (a n p) n under the assumption and combining flB.2j) we find 

(B.5) 
This result reads 



K\ Pc )^{a n p c T = e ^U\jl , * 



i 

n + 1 



a nPc 



2J1 



4A , 



23 



re— 1 1 



e 2 n 2 f 2r v 2 ( f 2 \ n+1 / f 2 \ " H 

K 7 n)^- (B.6) 



In 



2 2p 2 \ 2 J \ 4Ay 

The assumption (1B.3|) is satisfied if 



+ ^ » ^p- (B.7) 

Let us remember that a value of a" 1 is proportional to 2n for large n and v^A for large A and 
7„ can be expected to be 7„ ~ 1 although 7 n depends on A, e and n slightly. Sufficiently large 
n, large A and small e satisfy the above and therefore in several limits we find 

lim^»=-l, lim^ = -l, lim^^i = -l, (B.8) 

n— >oo 2 A— s-oo 2 e— >0 2 

whereas the case with n — 1 gives an important exception as 

lim— 1 = -1 + ai x 7 1 | e _ >0 . (B.9) 
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